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A statistical model for the parton distributions in the nucleon has proven its efficiency in 
the analysis of deep inelastic scattering data, so we propose to extend this approach to the 
description of unpolarized fragmentation functions for the octet baryons. The characteristics 
of the model are determined by using some data on the inclusive production of proton and 
' A in unpolarized deep inelastic scattering and a next-to-leading analysis of the available 

experimental data on the production of unpolarized octet baryons in e + e~ annihilation. 
Our results show that both parton distributions and fragmentation functions are compatible 
with the statistical approach, in terms of a few free parameters, whose interpretation will be 
discussed. 
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O . I. INTRODUCTION 

o, 

Following the first evidence of partonic substructure of the nucleon, by means of deep inelastic 
scattering (DIS), a large amount of experimental data have been collected in order to understand 
the parton structure of the nucleon. It is currently known that the Q 2 evolution behavior of the 
parton distribution functions (PDF) is well described by perturbative QCD. However, the parton 
Q H ' distributions at an initial scale reflect the nonperturbative quark and gluon dynamics of QCD bound 

states, which cannot be determined from first principles. For this reason, many parametrizations 
have been proposed jH H, 0, E( , but most of these expressions involve a large number of free 



parameters, with no clear physical meaning. However, some efforts have been recently made 
to parametrize unpolarized and polarized nucleon PDF based on a statistical approach 0, S|- The 
nucleon PDF, which involve in this framework a small number of free parameters, can well describe 
all available experimental unpolarized and polarized DIS data, so this is perhaps an indication that 
the PDF retain some important statistical features of the nucleon. The statistical approach of the 
nucleon PDF allows to make predictions, which were tested recently in a satisfactory way, by 
various experimental data in DIS 9j. In order to extend further this framework to new areas, it is 
natural to envisage, for example, its application to the description of the PDF of the other octet 
baryons. Unfortunately, these PDF are not directly accessible because, due to their short lifetimes, 
the hyperons cannot be used as a target in a DIS experiment. 

However, it is well known that the quark distributions in a hadron h, qh{x,Q 2 ) are related to 
the corresponding quark to hadron fragmentation function (FF) D^(x,Q 2 ), by means of the so 
called Gribov-Lipatov relations 0). This connection between two basic quantities of the hadron 
structure has been used in several recent works 11, 3], as an attempt to improve our present poor 



knowledge on the hadron FF. So this is our motivation to extend the statistical approach to a global 
description of the octet baryon FF and to check its validity against the available experimental data. 
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The paper is organized as follows. In Sect. 2, we review the main points of the framework and 
we give the basic procedure for the construction of the octet baryons FF in the statistical approach, 
in terms of Fermi-Dirac distributions. In Sect. 3, we determine all free parameters of the model 
by using some data on the inclusive production of proton and A in unpolarized deep inelastic 
scattering and a next-to-leading (NLO) fit to the available experimental data on the production 
of unpolarized octet baryons in e + e~ annihilation. The results of the fit for the cross sections and 
the obtained FF are also presented and we compare our approach with some previous works on 
baryon FF 0, Q, HI, Q . Finally in Sec. 4, we give our conclusions. 



II. THE STATISTICAL APPROACH FOR UNPOLARIZED FF 

In the statistical approach, a hadron can be viewed as a gas of massless partons (quarks, 
antiquarks, gluons) in equilibrium at a given temperature, in a finite size volume. In hadron 
production, when a parton fragments into a hadron, it picks up other partons from the QCD 
vacuum in order to form a specific hadron. The formation probability of the hadron is characterized 
by its statistical properties. Therefore we believe that statistical features which were proposed to 
build up the nucleon PDF in Ref. can be used also to construct the FF for the octet baryons. 
We assume that the parton (p) to hadron (h) FF Dp(x), at an input energy scale Qq, is proportional 
to 



[exp[(x - X )/x] ± l]- 1 , (1) 

where the plus sign for quarks and antiquarks, corresponds to a Fermi-Dirac distribution and the 
minus sign for gluons, corresponds to a Bose-Einstein distribution. Here Xq is a constant which 
plays the role of the thermodynamical potential of the quark hadronization into a hadron and x is 
the universal temperature, which is assumed to be equal for all octet baryons. It is reasonable to 
take its value to be the same as for the nucleon PDF, i.e. x = 0.099, according to Ref. @|. The 
statistical approach for the PDF allows to construct quark distributions of a given helicity and it is 
also possible to relate simply quark to antiquark distributions, resulting from chiral properties of 
QCD, as explained in Ref. p. All these physical quantities were determined precisely due to the 
existence of a hudge amount of data in unpolarized and polarized DIS. In the case of the octect 
baryon FF, the situation is different and the sarcacity of the polarized data does not allow such 
a clear separation, so we will restrict ourselves to the determination of the unpolarized quark FF, 
although the extention to the polarized case can be easily done. For the quarks q = u,s,d the FF 
are then expressed as 



A^Xirx" 



D (x, Q ) — - R , (2) 

y exp[(x — X^)/x\ + 1 

where X^ is the potential corresponding to the fragmentation q —* B and Qq is an initial scale, 
given below in Table 1. We will ignore the antiquark FF D^, which are considered to be strongly 
suppressed. The heavy quark FF Dq(x, Qq) for Q = c,b, i, which are expected to be large only in 
the small x region (x < 0.1 or so), are parametrized by a diffr active term with a vanishing potential 
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The initial scale Qq, which is flavor dependent in this case, is given below in Table 1 1 . This FF for 
Q —* B depends on b and a normalization constant A® for each baryon B. For the other quarks, 
we make some reasonable assumptions in order to reduce the number of parameters in addition to 
b, the universal power of x in Eq. (j2j). First we have the obvious constraints, namely, = D d 
for B = p, A. Moreover we assume that we need only four potentials, two for the proton X? = X v d 
and Xf and two for the hyperons = Xj and Xj where Y = A, X^, H~. Finally for the gluon 
to baryon FF Dg(x, Q 2 ), which is hard to determine precisely, we take a Bose-Einstein expression 
with a vanishing potential 

AB'b+l 

Df(x,Q 2 )= 7* (4) 
y exp(x/x) — 1 

We assume it has the same small x behavior as the heavy quarks and it is the same for all baryons. 
The normalization constants A^ : Ag and Aq will have also to be determined by fitting the data, 
a procedure we present now. 

III. DETERMINATION OF THE PARAMETERS FROM THE DATA ANALYSIS 

The experimental data which was used is twofold. First,Jiadron production in DIS gives access 



to a direct determination of DP(x, Q ) [ljfl and D£(x, Q ) |18(], in a limited range of x and Q . We 



have also determined the free parameters of the model by making an analysis of the differential 
cross section for the semi-inclusive hadron production process e + e~ — > h+X, for p, A, S, H. Here we 
note that in these experiments, they usually do not distinguish between B and B, so we will include 
both contributions in our calculations by making the natural assumption Dg(x, Q 2 ) = D?(x, Q 2 ). 
The differential cross section can be expressed as 

E ED, EH, 



i 

1 da _ ^ f drj ^hfXE 2 \ 1 



Vtot dXE 



„ J V i] atot 



where xe = 2Eh/y/s, which is the energy Eh of the produced hadron scaled to the beam energy 
Q/2 = y/s/2. 2 The subscripts T and L denote the contributions due to transverse and longitudinal 
polarizations, respectively. The summation q includes quarks, antiquarks and gluon contributions. 
In the following we will set the renormalization and the factorization scales equal to [i 2 R = fsp = Q 2 . 
The definitions of the functions C\ and C\ follow those of Ref. 21]. In Eq. (fKj). a^t is the total 
cross section for the process. 

Aira 2 , /, a As 



^ = ^£^(^l + ^j, (6) 

where N c is the color number and we shall take N c = 3, a is the QED fine structure constant and 
a s (s) is the strong coupling constant. The sum in the above equation should be over all active 
quarks and antiquarks. The electroweak charges in Eq. © can be expressed as 

<5g0) = + 2v e v q e q pi{s) + (a 2 e + v 2 )(a 2 + v 2 )p 2 (s) , (7) 



1 Due to the fact that the input scale of the t quark is above the highest energy data investigated in this work, it 
does not contribute to our analysis. 

2 When experimental data are presented in the momentum scaling variable z p = 2ph/Q, we made a conversion to 
xe in the following figures. 
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with 



1 s(M\ - a) . . 

P1[S} ~ 4 sin 2 6 W cos 2 % (Mf - s )2 + Mfr| ' ( j 



1 s 2 
P2{S) = 16 sin 4 9 W cos 4 W (Mf - s )2 + Mfr| ' (9) 

o e = -l, a g = T 3g , (10) 

f e = + 2sin 2 (9 1 y, u g = T 3q - 2e q sin 2 W , (H) 

where e q is the electric charge of the quarks and (vj), (a^) the electroweak vector and axial couplings 
of the electron and the quarks to the Z, respectively. T 3q is the third component of the strong 
isospin, Tz q = 1/2,-1/2 for up-type quark and down-type quark, respectively. 6\y is the weak- 
mixing angle, and Mz, ^ z are the mass and width of the Z. 

The FF have been evolved at NLO following a method defined in Ref. [f|, where for the DGLAP 
equations we used the timelike splitting functions calculated in Refs. j^ . For baryon production in 
e + e~ collisions, we made a NLO fit with formula © of the experimental data from Refs. 24|-[42T|. 



where we restricted xe > 0.1. In addition, the input scales choice |2l| and A(MS) are given in 
Table HI 

TABLE I: Input scale Q and A(MS) in GeV unit. 



quark 


u,d,s 


c 


b 


t 


Qo 


0.632 


1.4 


4.5 


175 


A(MS) 


0.299 


0.246 


0.168 


0.068 



Now, let us report the values of the free parameters we have obtained from the NLO fit: 

XI = 0.648, = 0.247, X A = 0.296, X S A = 0.476 . . 

b = 0.200, b = -0.472, = 0.051. 

These parameters have similar values to those obtained for the nucleon PDF in particular, for 
the thermodynamical potentials, so the intrinsic properties of the quarks when observed in DIS 
or in fragmentation processes seem to be preserved. Notice that in the nucleon PDF the u quark 
which is dominant has the larger potential and here we have analogously, X? > X P and X A > X A , 
a situation which is natural to expect. The other parameters A qi , A q2 , for the quark to baryon 
FF, and Aq for the heavy quarks are given in Table ITU together with the quark content q±, q2 
and Q for each baryon B. Clearly these normalization constants are decreasing when going from 
the proton to the heavier hyperons, following the magnitude of the corresponding measured cross 
sections. 

Finally, let us comment on our results. With the data set on proton, A, S + , £~ and H~ 
production, we get a x 2 = 227.5 for 206 experimental points (x% > 0.1), i.e. xVpoint = 1.1 for 
our NLO fit. This agreement of our results with experimental data is satisfactory and it confirms 
that the statistical approach is also successful in the description of the octet baryons FF. 
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TABLE II: Values of the normalization constants of the FF for the octet baryons 





Baryon 


9i 


92 


A B 




A B 
A Q 


p(uud) 


u = d 


s 


0.264 


1.168 


2.943 


A(uds) 


u = d 


s 


0.428 


1.094 


0.720 


T, + (uus) 


u 


s 


0.033 


0.462 


0.180 


Y.-(dds) 


d 


s 


0.030 


0.319 


0.180 


S~(dss) 


d 


s 


0.023 


0.082 


0.072 



In Figs. n& n d[21we display the available DIS data, which yield directly DP(x, Q 2 ) and D^(x, Q 2 ), 
and the result of our fit. In Figs. 13161 we show a comparison of the calculated cross sections for 
various baryons with experimental data. In Fig. |31 we give our results for the proton cross section 
in electron-positron annihilation at y/s = 22,29,34 GeV and at the Z-pole i.e. y/s = 91.2 GeV. 
The agreement with data is very good even in the low xe region {xe < 0.1) as shown on the same 
figure, although these data points were not included in our fit. 

For the A production, our results are displayed in Fig. ^ for the energies y/s = 
14,22,29,33.3,34.8,42.1,91.2 GeV. It is clear from Figs. El and H which display a sizeable en- 
ergy domain, that the scaling violations in this range are very small. In order to improve the small 
xe behavior, we need to modify the evolution of the FF in this region, and also include finite 
mass corrections and modifications of the splitting functions (see a discussion in Ref. but 
such corrections are outside the scope of the present paper. In Fig. GD we give our results for E 1 * 1 
production at y/s = 91.2 GeV and in Fig.|f)]our results for H~ production at the Z-pole. For these 
strangeness production processes, our results show a satisfactory agreement. The corresponding 
quark to baryon FF are presented in Fig. [7| For all hyperons the strange quark FF dominates 
largely over the u, d quarks, which seems natural. For A, we are in agreement with a model for 
SU(3) symmetry breaking in Ref. which leads to ~ 0.07-D^. In Ref. one also finds 
smaller than and Dq is strongly suppressed. This is in contrast with the situation of Ref. 
[13j, where u, d and s are assumed to contribute equally. However the heavy quarks have a pattern 
similar to Ref. ^3], with a sizeable contribution only for x < 0.1 and a fast dropping off for large 
x. This is also at variance with Ref. [ifj], where D^/D^ decreases from 1 to 0.2 when x goes 
from to 1. For the proton it is surprizing to see that the u-quark FF dominates only at large x, 
whereas the strange and heavy quarks contribute substantially for x < 0.3 or so. Finally we notice 
that for the A, the heavy quarks have a pattern similar to Ref. [l3T |. with a sizeable contribution 
only for x < 0.1 and a fast dropping off for large x. 

IV. CONCLUSIONS 

With the motivation to check whether the octet baryons FF have similar statistical features as 
the nucleon PDF previously studied, we extended the statistical approach to analyze some data on 
the inclusive production of the octet baryons. We found that these FF can be well described with 
a small number of free parameters, whose interpretion was discussed. We obtained a satisfactory 
description of the unpolarized experimental data, suggesting that the statistical approach, with 
Fermi-Dirac type FF, works equally well, compared to other parametrization forms. The semi- 
inclusive DIS data for proton and A production, which give strong constraints, have allowed a 
flavor separation between u, d and s quarks FF, an interesting result which remains to be more 
seriously checked in the future. 

In our present analysis, we did not introduce polarized FF, although our formalism can be 
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easily extended to this case. Actually, it is impossible to extract some reliable information on the 
polarized FF for the moment, due to the scarcity of the data, but we hope this will be possible in 
the future. 
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FIG. 2: The FF for u quark to A, D^x.Q 2 ), as a function of x at Q 2 — 2.5GeV 2 . The experimental data 
are from Ref. [Isj . 
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FIG. 3: Cross sections for proton production in e + e annihilation at several energies as function of Xe- The 
experimental data are from Refs. jil El El ElES El] . 
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FIG. 4: Cross sections for A production in e + e annihilation at several energies, as function of xe- The 
experimental data are from Refs. E3, EH EE EE EE EE E3, EE Ii3 
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FIG. 5: Cross sections for E ± production in e + e annihilation at the Z-pole as function of Xe- The 
experimental data are from Ref. |4lj . 
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FIG. 7: The quark to octet baryons FF D^(x,Q 2 ) and D^{x,Q 2 ) (B = p, A, E ± , E~ , q = u,d,s and 
Q = c, b, t), as a function of x at Q = 91.2GcV. Note that we used different vertical scales in the upper and 
lower parts of the figure. 



